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ABSTRACT 

We revisit the concept of a blazar sequence that relates the synchrotron peak frequency (/'peak) in 
blazars with synchrotron peak luminosity (L peak , in vh v ) using a large sample of radio-loud AGN. 
We present observational evidence that the blazar sequence is formed from two populations in the 
synchrotron /'peak — L pea k plane, each forming an upper edge to an envelope of progressively misaligned 
blazars, and connecting to an adjacent group of radio galaxies having jets viewed at much larger 
angles to the line of sight. When binned by jet kinetic power (Lk in ; as measured through a scaling 
relationship with extended radio power), we find that radio core dominance decreases with decreasing 
synchrotron L pea k, revealing that sources in the envelope are generally more misaligned. We find 
population-based evidence of velocity gradients in jets at low kinetic powers (~ 10 42 — 10 44 5 erg s" 1 ), 
corresponding to FR I radio galaxies and most BL Lacs. These low jet power ‘weak jet 7 sources, 
thought to exhibit radiatively inefficient accretion, are distinguished from the population of non- 
decelerating, low synchrotron-peaking (LSP) blazars and FR II radio galaxies (‘strong 7 jets) which 
are thought to exhibit radiatively efficient accretion. The two-population interpretation explains the 
apparent contradiction of the existence of highly core-dominated, low-power blazars at both low and 
high synchrotron peak frequencies, and further implies that most intermediate synchrotron peak (ISP) 
sources are not intermediate in intrinsic jet power between LSP and high synchrotron-peaking (HSP) 
sources, but are more misaligned versions of HSP sources with similar jet powers. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: jets 


1. INTRODUCTION 

The success of the current model for radio-loud Ac- 
tive Galactic Nuclei (AGN) as accreting super-massive 
black holes with bipolar, relativistic jets has rested on 
the ability of this physical paradigm to successfully ex- 
plain the chameleon nature of the spectral energy dis- 
tribution (SED) of this family of extreme objects. The 
jet model is inherently anisotropic and allows us to unify 
the highly luminous and variable blazars as relativisti- 
cally beamed counterparts of Fanaroff-Riley (FR) radio 
galaxies (Fanaroff & Riley 1974). 

The typical blazar spectrum exhibits a double peak 
(in uF v ), attributed to synchrotron and inverse Compton 
(IC) radiation from the relativistic jet, which dominates 
in blazars due to the alignment of the jet axis near to our 
line-of-sight (for a review of leptonic blazar models, see 
Bottcher 2007). The synchrotron component peaks from 
sub-infrared (IR) energies up to hard X-rays, with the IC 
component peaking at much higher energies (typically 
gamma-rays). Spectroscopically, blazars are tradition- 
ally divided between optically ‘featureless’ BLLac ob- 
jects and (broad) emission line-emitting Flat Spectrum 
Radio Quasars (FSRQ). The long-standing orientation- 
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based unification scheme for radio-loud AGN proposes 
BL Lacs and FSRQ as counterparts to FR I and FR 
II radio galaxies, respectively, based on similar spectra, 
morphologies, and range in extended radio luminosity 
(Urry & Padovani 1995). 

However, violations to this scheme are well-known, 
and findings of powerful, FR H-like BL Lacs, and 
low-power FSRQ (Landt et al. 2006; Landt & Bignall 
2008; Kharb et al. 2010) appear to breaJk the sim- 
plest dichotomy (going back to the original FR classi- 
fication) between the lobe-dominated, high-power FR 
II galaxies and the low-power, edge-darkened FR Is. 
Kharb et al. (2010) also found BL Lacs displaying SEDs 
and hot-spots typical of quasars while others have ob- 
served BL Lacs exhibiting broad lines in low contin- 
uum states (Vermeulen et al. 1995; Corbett et al. 2000; 
Raiteri et al. 2007) . 

Further, the importance of optical spectral type in 
connecting radio galaxies with their aligned counter- 
parts is difficult to assess due to the strong effects 
of relativistic beaming on equivalent width measure- 
ments on which the blazar classifications are based 
(Georganopoulos & Marscher 1998; Ghisellini et al. 
2010). At the same time, the designation of Low and 
High Excitation Radio Galaxies (LERG and HERG, cor- 
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responding to prominence of emission lines as measured 
with an excitation index, see Hine & Longair 1979) is 
increasingly found to be mixed among FR I and FR II 
morphological types; all of this creates serious problems 
for a unification scheme which connects spectral type 
with the morphology. 

1.1. The Blazar Sequence 

The synchrotron peak frequency (^ pea k) may take on 
a wide range of values (from 10 12 to 10 18 Hz at the ex- 
tremes), and is one of the principle ways to classify indi- 
vidual blazars. While FSRQ sources are typically found 
to have low i/ pea k (< 10 14,5 Hz), BL Lacs are found to 
span the entire range, with low, intermediate, and high- 
peaking BL Lacs (LBL, IBL, and HBL) typically dom- 
inating surveys selecting in the radio, optical, or X-ray, 
respectively. We adopt the generic terms for low, inter- 
mediate, and high synchrotron-peaking (LSP, ISP, HSP) 
blazars independently of the spectroscopic type, as in 
Abdo et al. (2010c). 

Fossati et al. (1998) found an anti-correlation between 
the frequency of the synchrotron peak and the bolomet- 
ric and synchrotron peak luminosities and suggested a 
link between the power of the source and 2/ pea k. This 
sequence appeared to match a spectroscopic one, from 
predominantly powerful FSRQ sources at the low-peak, 
high-luminosity end through LBLs and finally HBLs, at 
the low-luminosity end. Ghisellini et al. (1998) suggested 
that more efficient cooling of particles in the jets of high 
luminosity blazars is responsible for the lower peak fre- 
quencies. However, it is not clear how the continuous 
blazar sequence fits in with the bi-modality of morphol- 
ogy and spectral type in radio galaxies. 

1.2. Beaming and Quantitative Unification 

Orientation-based unification predicts that misaligned 
blazars will drop in luminosity and synchrotron peak 
frequency according to the decrease in Doppler boost- 
ing with increasing angle. In the simplest case (a single 
Lorentz factor and a convex spectrum), the ratio of de- 
crease in log frequency to log luminosity will be approx- 
imately 1:4, roughly orthogonal to the sequence. (This 
holds for moving features in the jet; for a standing shock 
model the relation is 1:3, see e.g., Lind & Blandford 
1985). The connection between blazars and the nuclear 
emission in radio galaxies can naturally inform our jet 
models. The synchrotron nature of the nuclear radio 
galaxy emission has been confirmed by comparison of 
the broad-band spectral indices and SED characteris- 
tics of FR I nuclei with BL Lacs (Chiaberge et al. 1999; 
Capetti et al. 2000; Chiaberge et al. 2000; Capetti et al. 
2002; Trussoni et al. 2003), where the effects of relativis- 
tic beaming (rather than obscuration) accounted for the 
10— 10 4 difference in radio and optical luminosities be- 
tween the BL Lacs and radio galaxies. 

Capetti et al. (2002) and Trussoni et al. (2003) agree 
upon jet bulk Lorentz factors of T ~ 4-6 in order to 
account for the difference in radio galaxy and blazar 
luminosities, while Hardcastle et al. (2003) found T ~ 
3 sufficient to match the luminosity functions of the 
B2 and EMSS and 1 Jy samples. These values agree 
with VLBI studies of the pc-scale jets of HSP BL Lacs 
which show that they are barely super-luminal, with av- 
erage bulk Lorentz factors T ~ 2-3 (Piner et al. 2001; 


Edwards &: Piner 2002; Giroletti et al. 2004; Piner et al. 
2008, 2010); however, several LSP BL Lacs show high ve- 
locities on VLBI scales (Jorstad et al. 2001; Cohen et al. 
2007). Interestingly, Hardcastle et al. (2003) also found 
that they could not reproduce the redshift distribution 
of 1 Jy BL Lac sample unless some of the high-redshift 
BL Lacs corresponded to FR IIs. 

In contrast, high Doppler factors (~ 50 or greater) 
must be invoked to explain the brightness and rapid 
variability of TeV blazars, as the jet must be highly 
relativistic in order to avoid absorption by co-spatial 
IR photons (Dondi & Ghisellini 1995; Aharonian et al. 
2007; Ghisellini et al. 2009; Giannios et al. 2009). The 
discrepancy in apparent jet speeds and radio galaxy nu- 
clear emission can be resolved at least for FR I/BL 
Lac sources by invoking the presence of structured 
jets with multiple emitting regions of varying speeds. 
This includes the spine-sheath axisymmetric model ex- 
plored by Swain et al. (1998), Chiaberge et al. (2000), 
and Ghisellini et al. (2005), and the decelerating model 
of Georganopoulos & Kazanas (2003), all of which imply 
a high-speed spine or base which dominates the flux when 
the source is aligned, and a slower component which dom- 
inates if the source is misaligned. Other blazar stud- 
ies also support the decelerating model for low power 
blazars, as recent VLBI observations of the MOJAVE 
sample 1 reveal that the these sources also have lower ap- 
parent jet speeds, versus the high-power FSRQ which 
maintain collimated jets and consistent jet speeds over 
longer distances (Lister et al. 2009). 

1.3. The Blazar Envelope 

The simplest scenario behind the appearance of a 
blazar sequence is that the physical properties of extra- 
galactic jets (as measured through the synchrotron peak 
frequency and luminosity) are a function of a single pa- 
rameter, the jet kinetic power (Ljdn)- If there is a one-to- 
one correspondence between ^ pea k and L^n, the original 
blazar sequence could be seen as the extreme ‘edge’ of 
an envelope, consisting of the most aligned (and there- 
fore brightest) objects. More misaligned sources should 
appear below the sequence, causing it to widen. If the 
jets are monoparametric in terms of jet power, these 
misaligned sources should fall in a line away from their 
aligned counterparts on the sequence, according to the 
drop in Doppler factor, forming ‘tracks’ of similar kinetic 
power in the envelope of blazars populating the ^ pea k - 
L p eak plane. 

The original blazar sequence was based on 126 objects, 
a few dozen of which were detected in gamma-ravs. Over 
a decade later there are thousands more blazars and 
blazar candidates in the literature and a much greater 
density of multi-wavelength data in general. Many of the 
newer samples have been designed to combat the early se- 
lection effects of the radio and X-ray surveys, such as the 
low radio power CLASS sample (Caccianiga & Marcha 
2004), and it is now clear that as our surveys go deeper, 
many more low-luminosity blazars are emerging. These 
newer samples challenge both the sequence and its sim- 
plest extension that considers the sources in the enve- 

1 MOJAVE (Monitoring Of Jets in Active galactic nuclei with 
VLB A Experiments): Over 250 bright blazars being continuously 
monitored with VLBI. (www.physics.purdue.edu/mojave) 
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TABLE 1 
Blazar Samples 


Sample 

Abbrev. 

No. of 
Objects^ 

References 


1 Jansky blazar sample 

1 Jy 

34 

Stickel et al. (1991) 


2-iegree field (2dF) QSO survey 

2QZ 

56 

Londish et al. (2002, 2007) 


2 Jansky survey of flat-spectrum sources 

2 Jy 

213 

Wall ic Peacock (1985) 


The Candidate Gamma-Ray Blazar Survey 

CGRaBS 

1625 

Healey et al. (2008) 


Cosmic Lens All-sky Survey 

CLASS 

236 

Marcha et al. (2001); Caccianiga et al. 
Caccianiga & Marcha (2004) 

(2002); 

Deep X-ray Radio Blazar Survey 

DXRBS 

283 

Perlman et al. (1998); Landt et al. 

Padovani et al. (2007) 

(2001); 

EGRET blazar sample 


555 

Mukherjee et al. (1997) 


E:nstein Slew survey sample of BL Lac objects 


48 

Perlman et al. (1996) 


Fermi blazar sample 

1FGL 

689 

Abdo et al. (2010a) 


Hamburg-RASS bright X-ray AGN sample 

HRX 

104 

Beckmann et al. (2003) 


Parkes quarter- Jansky flat-spectrum sample 


878 

Jackson et al. (2002); Hook et al. (2003); Wall et al. 

Radio-Emitting X-ray source survey 

REX 

55 

Caccianiga et al. (1999) 


Radio-Optical-X-ray catalog 

ROXA 

816 

Turriziani et al. (2007) 


RASS - Green Bank BL Lac sample 

RGB 

127 

Laurent-Muehleisen et al. (1999) 


E'instein Medium-Sensitivity Survey of BL Lacs 

EMSS 

44 

Morris et al. (1991); Rector et al. (2000) 


RASS - SDSS flat-spectrum sample 


501 

Plotkin et al. (2008) 


Sedentary survey of high-peak BL Lacs 


150 

Giommi et al. (1999, 2005); Piranomonte et al. (2007) 


* Number of sources listed is total sources in sample (not necessarily number of blazars). 


tended radio luminosity are related. We therefore chose 
the low-frequency extended luminosity (also commonly 
referred to as steep or lobe emission) at 300 MHz (L 300 ) 
as an estimator of the jet power for this work. This value 
is then converted to Lkin from the linear fit (with l-a er- 
ror), 

log Lkin = 0.64(±0.09) (log L300 - 40) 

+ 43.54(±0.12) [ergs s -1 ] 

taken from Cavagnolo et al. (2010). Although the scatter 
in the scaling is not insignificant (~ 0.7 decade in Lkin), it 
is adequate for our purpose of classifying sources in broad 
bins of jet power. This value (Lkin) is what we generally 
refer to throughout the paper as the jet (kinetic) power, 
not to be confused with L pea k or other observables of the 
jet which are affected by beaming. 

2 . 2 . 1 . Extended Radio Emission 

The extended radio emission from the lobes of radio 
galaxies differs from the beamed core emission in be- 
ing both non-variable and generally exhibiting a much 
steeper spectrum. Thus one way to isolate the extended 
emission is to look for a change in the spectral index 
at low frequencies. For our sample, a spectral decom- 
position (SD) method was applied to all sources which 
had data points at five or more frequencies below 10 11 
Hz. Sources were tested for a two-component spectrum 
by fitting alternately with a double power-law and then 
a single power law (assuming normal error). Because 
the latter can be expressed as a particular case of the 
former, this allowed us to use the maximum likelihood 
(L) ratio test for nested models. In this test, the statis- 
tic — 2Alog(L) distribution can be approximated as a x 2 
with ri 2 - n\ degrees of freedom (Wilks 1938). The double 
power-law model has 712 = 4 degrees of freedom (slope 
and normalization for each power law), while the null hy- 
pothesis model has n\ = 2, for a total of two degrees of 


freedom. We rejected the two power-law fit for a signif- 
icance level a >0.05. The best steep-spectrum fit from 
the hypothesis test was then used to calculate L 300 - 

The normal approximation for the spread around the 
power law is a rough approximation, and for eight sources 
it appears the fit failed to find a significant two-powerlaw 
solution due to a great deal of variability in the radio 
core points. L 300 was estimated for these sources by 
fitting hand-selected low-frequency radio points corre- 
sponding to a stable emission. component persistent be- 
neath the core variability. This method was also used 
for a few sources which passed the test, where core vari- 
ability caused the steep power law fit to be somewhat 
more flat than would be fit by eye (a correction of a few 
percent). We note that inclusion of these eight sources 
and the corrections do not materially affect the results 
presented. 

The SD method is the most reliable way to separate the 
boosted core emission from the isotropic lobe emission, 
as it is model independent. Another method is to mea- 
sure the extended emission from core-subtracted VLA 
images (typically at 1.4 GHz). In order to expand our 
sample and counter some selection bias (see longer dis- 
cussion in §4.2), extended flux measurements were taken 
from the literature where available (see Table 2 for ref- 
erences). One difficulty with these measurements is that 
an assumption must be made about the spectral index 
used to extrapolate from 1.4 GHz to the 300 MHz lumi- 
nosity. It is also clear from sources with both a visible 
steep spectrum and VLA data that these are not always 
in agreement. Even for the same object, VLA estimates 
by different authors can vary widely, as in the case of 3C 
446, with 1.4 GHz extended flux estimates of 91.6, 271, 
and 4030 mjy (Kharb et al. 2010; Cooper et al. 2007; 
Antonucci &; Ulvestad 1985), none of which agrees with 
a clear steep component which suggests a flux ~ 800 mjy 
at that frequency. However, examination of the over 100 
sources with both VLA and SD estimates reveals that 
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lope simply to be de- beamed blazars. Contrary to what 
is anticipated by the sequence, Landt &: Bignall (2008) 
showed that HSP and LSP BL Lacs objects cover similar 
ranges in extended radio power. Caccianiga & Marcha 
(2004) found several low luminosity, low i/ pea k sources 
with relatively high core dominance compared to brighter 
objects at similar low y pea ki suggesting that the former 
are not simply related to the latter through misalign- 
ment, as one might naively expect. 

In this paper we examine the available evidence for 
a ‘blazar envelope’ beneath the upper-limit of the most 
aligned blazars, due primarily to mis-alignment. We test 
the possibility that jets are monoparametric engines by 
examining how the source power is related to other prop- 
erties of the blazar such as location in the v peak - L pea k 
plane, and re-evaluate the blazar sequence in light of 
the great increases in multi-wavelength, multi-epoch data 
available for blazar and radio galaxy jets. 

The paper is organized as follows: in §2 we explain 
the sample selection and describe our statistical and 
phenomenological models which allowed us to estimate 
the important envelope parameters for blazars and radio 
galaxies. In §3 we present the results of the first round 
of data collection and analysis and the ‘blazar envelope’. 
In §4 we discuss the implications of the blazar envelope 
and the potential for constraining jet models. In §5, we 
summarize our findings. Throughout the paper we adopt 
a cosmology based on the concordance model with Qm 
= 0.27, Ha = 0.73, and Hq = 71 km s -1 Mpc” 1 . Energy 
spectral indices are defined such that f u oc v~ a : and all 
luminosities (L u ) are references to powers (i/L„ implied). 

2. METHODS 

In order to investigate the nature of the blazar enve- 
lope, v:e require a sample of blazars with a well-sampled 
SED such that a reliable estimate of the synchrotron peak 
location can be made. The scarcity of such blazars may 
recommend utilizing as many sources as are known with- 
out regard for flux-limited samples, but this approach 
risks losing the ability to understand the limits of our ob- 
servations. We compromise by combining multiple well- 
defined (flux-limited) samples of flat-spectrum sources. 
The overall sample definition is thus complex, but does 
allow for the examination of biases by population mod- 
eling, as discussed in §4.2. 

2.1. The Blazar Sample 

The flat spectrum candidate sample was made from the 
surveys listed in Table 1. A multi- frequency database of 
all sources in this master list was compiled by a combina- 
tion of all the data in the original publications, with ad- 
ditional data from the literature and from the SEVIBAD 
and NED online databases, for a total of 3773 sources 
(taking into account duplicates). For 69 sources, previ- 
ously unpublished Swift data were available which were 
utilized to produce fluxes and X-ray indices. Spectro- 
scopic classification was taken from NED, SIMBAD, or 
the most recent publication, and non-blazar sources were 
removed, leaving 3014 blazars. A further reduction was 
made by removing sources with sparse multi-frequency 
SED sampling. The minimum required sampling was 
photometry at 3, 2, and 1 or more unique frequencies for 
the (rest frame) frequency ranges < 10 11 5 , 10 11,5 -10 16 , 


and > 10 16 Hz, respectively, yielding 1327 sources suit- 
able for further analysis. We discuss the properties of 
the discarded sources in §4.2.2. 

For sources with unknown redshift (~ 16% of the sam- 
ple), we use the relationship between optical host galaxy 
magnitude and redshift of Sbarufatti et al. (2005) to es- 
timate a value for these sources, taking the lower limit 
when exact host magnitudes were unavailable. All lu- 
minosities and frequencies were corrected to rest-frame 
values before spectral fitting; the global fits (rather than 
an assumed index) were then used to find k-corrected 
rest-frame quantities of interest. 

2.2. Jet Power 

Measuring the jet kinetic jet power of our blazars and 
radio galaxies allows us to see if the areas occupied by jets 
of similar power in the i^ pe ak — L pea k plane can be under- 
stood simply as an orientation effect (as it would be an- 
ticipated in the original blazar sequence), or if additional 
considerations are needed. We note that for our purposes 
we need only classify sources into a few broad bins of jet 
power (as we show below, a width of one decade gives 
us a sufficient number of well populated jet power bins), 
and an extremely accurate estimation of the jet kinetic 
power is not needed. 

Rawlings & Saunders (1991) estimated jet power by di- 
viding the lobe energy content of radio galaxies by their 
age. The energy content was estimated through mini- 
mum energy arguments ( e.g . Miley 1980) derived from 
measurements of the low frequency radio emission where 
the lobe emission is assumed to dominate over the jet, 
while the source age was estimated from spectral breaks 
in the lobe spectra through synchrotron aging arguments. 
Willott et al. (1999) revisited the issue of obtaining the 
jet power from the low frequency extended luminosity. 
They studied in detail the parameters that come into this 
relation, including the unknown energy content of the 
non-radiating particles and the time evolution of the lobe 
luminosity (see also Blundell et al. 1999) and presented 
a scaling relation which included a careful parametriza- 
tion of the various uncertainties (see their Section 4.1, 
and Equation 11). 

The actual normalization for this scaling between the 
jet kinetic power and the low frequency extended lumi- 
nosity (L ex t) came from a very different route. The dark 
cavities seen in {e.g., Chandra) X-ray images are created 
by AGN jet activity and can be used to estimate the pAV 
work needed to inflate them. Importantly, this energy 
estimate does not depend on the minimum energy as- 
sumption and does not require knowledge of the amount 
of non-radiating matter in the lobes. The jet power Lkin 
is estimated by dividing the pAV work needed to ex- 
pand the cavity by the cavity age found by estimating the 
lobe buoyancy time (Birzan et al. 2008; Cavagnolo et al. 
2010 ). 

While this method is limited to a small number of 
sources at present, the Lkin — L ext relation covers al- 
most six orders of magnitude in jet power, including 
both FR I and FR II sources. The slope of the scaling 
Cavagnolo et al. (2010) find is compatible within errors 
with that found by Willott et al. (1999) and requires a 
lobe energy content larger by at least a factor of 100 com- 
pared to that inferred by minimum energy arguments. 
The cavity work shows convincingly that Lki ;i and the ex- 
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Fig. 1. — SEDs for three sources in which the extended radio 
emission is clearly visible as a steep component at low frequen- 
cies.. The FSRQ QSO B0923+392 (top) has a very clear steep 
component at low frequencies, dominating up to ~ 10 9,5 Hz, and a 
prominent ‘big blue bump’ from the accretion disk, suggesting that 
it may be relatively misaligned compared with the Fermi-detected 
BL Lac IJy 2007+77 (middle), because they differ by ~ 2 decades 
in extended luminosity but only ^ 1 decade in beamed power. This 
latter source has a similar extended radio spectrum to the radio 
galaxy 3C 338 (bottom) , however in the radio galaxy more of the 
extended emission is visible due to the nuclear emission being much 
less beamed. We also show estimates of the extended emission at 
1.4 GHz as measured from VLA maps (diamonds), which generally 
agree with the spectral fits. The dashed curves are the best-fit of 
the blazar model described in the text and in Appendix A. 

typical differences are less than 0.5 in log luminosity, with 
the best agreement found using a low-frequency index of 
ol = 1.2 for the extrapolation to 300 MHz. 

The full SED for three sources is shown in Figure 1: the 
bright quasar QSO B0923+392 (top), the LSP BL Lac 
IJy 2007+77 (middle), and the FR I radio galaxy 3C 
338 (bottom). In all cases, extended radio emission was 
detected in our two-component SD test - these points 
are shown in lighter colors of pink, green, and orange, 
respectively. (Estimates of the extended luminosity de- 
rived from VLA maps are also shown as diamonds). We 
note that while IJy 2007+77 and 3C 338 have similar ex- 
tended luminosities at 300 MHz, the much greater core- 
dominance of the blazar is evident, and results in much 
less of the extended emission being visible. Also, while 
the quasar is evidently more powerful than the other two 
sources by two orders of magnitude in L 3 oo, we see that 
it is likely more misaligned than IJy 2007+77, from the 
greater visibility of the steep component and the promi- 
nent accretion disk bump in the optical. This relatively 
bright source was also not a significant detection in the 
Fermz/LAT 1-year catalog, unlike IJy 2007+77. This 
could be explained if the gamma-ray emission is more 
beamed than that at lower frequencies. 

2.2.2. Connection Between Extended and Unbeamed Core 
Emission 

For the SD method, it is important to note that failure 
to find a two power-law solution for a source can just as 
easily be due to a lack of low-frequency data coverage as 
to the source emission being highly jet dominated down 
to low frequencies. The approximate frequency at which 
the beamed emission becomes dominant varies from 10 6 
to 10 13 Hz in our sample. In Figure 2 we plot the ‘cross- 
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Fig. 2. — The crossing frequency, Across is the approximate fre- 
quency at which the SED becomes dominated by the beamed jet 
emission rather than the steep extended emission. Radio galaxies 
are shown as empty squares, and blazars as filled circles. As log 
core dominance (Rce) increases due to beaming this crossing point 
moves to lower energies. Many sources will remain core-dominated 
to low frequencies, even to the limit of ground-based observations 
at ~ 10 MHz. The broken line shown is the linear correlation be- 
tween the plotted variables (r = 0.87). The solid red line marks the 
300 MHz frequency at which the extended emission is measured. 

mg’ frequency z/ croas of the fitted jet SED and steep com- 
ponent and the log core dominance Rce (we define Rce 
as the log of the ratio of core to extended luminosity at 
1.4 GHz). The correlation is well- fit by a linear regression 
model (correlation coefficient r = 0.87). This correlation 
can arise from two coexisting causes: an intrinsic spread 
in unbeamed core power for a given extended luminosity 
and the effect of beaming. If the effect of the natural 
spread of core power dominates over the effects of beam- 
ing (as would be the case if beaming of the core emission 
is weak or absent), then we expect to see radio galaxies 
and blazars mixed along this correlation (this is because 
the range by which each particular source moves before 
being classified as a radio galaxy will be smaller than the 
initial intrinsic range that the aligned blazars populate). 
This is, however, not what we find, as shown in Figure 
2. Instead, we observe a clear separation of radio galax- 
ies and blazars which suggests that the effect of beaming 
dominates over the intrinsic spread and that the range of 
unbeamed core emission for a given extended luminosity 
be narrow. 


2.2.3. Caveats 

Many of our sources which lack VLA imaging have ob- 
servations down to very low frequencies (30 — 200 MHz), 
and show no signs of a steep component. It is impossible 
to tell if these sources are highly core-dominated to the 
lowest observed frequencies or if the lobe emission has 
a flat spectrum which is indistinguishable from the core 
emission. The latter scenario may be a concern for the 
very low power sources, as it is possible that the break 
frequency (Le., the point at which the lobe emission be- 
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Fig. 3. — Examples of typical blazar SEDs in the master sam- 
ple of flat-spectrum sources . Visual inspection of over 1300 SEDs 
was used to separate sources with unambiguous synchrotron peaks 
(top, 653 sources total) from those with under-sampled SEDs which 
might be consistent with several possible peak locations (bottom). 
The dark line is the best-fit model (described in Appendix A) for 
each case, while the dashed lines show other possible SED shapes 
consistent with the more sparsely-sampled SED. 

comes steep) may increase with decreasing synchrotron 
power as Brunetti et al. (2003) found for hot spots in ra- 
dio galaxies. For these sources it may be advisable to 
perform VLA map subtractions at multiple frequencies 
to confirm such a spectrum. 

Finally, we must point out the risk in the common 
practice of assuming that flux below ~ 300 MHz, taken 
without examination of the spectral characteristics, is 
actually from the extended isotropic component, rather 
than the core. This assumption would be clearly inap- 
propriate for many of our sources, as shown by Figure 2, 
and could lead both to an overestimate of the jet power 
and an underestimate of the radio core dominance often 
used to infer the degree of mis-orientation of the source. 

2.3. Fitting the Synchrotron Spectrum 

For the 1327 blazars making the initial cut based on 
SED sampling, we utilized a phenomenological paramet- 
ric model in order to estimate the synchrotron peak fre- 
quency and luminosity, The basic shape is parabolic, 
with a power-law extension to low frequencies and a ris- 
ing power-law to simulate the beginning of the IC spec- 
trum in the X-rays (Fossati et al. 1997). These com- 
ponents were free to vary within a reasonable range in 
shape or slope, respectively, and in relative normaliza- 
tion, for a total of 6 parameters. In sources where a 
big blue bump was suggested on visual inspection, this 
feature was added to the SED fit by adding a template 
accretion-disk spectrum (see Frank et al. 2002), adding 
one additional parameter for the bump peak luminosity 
(see Appendix A for a summary of the model). 

To fit our model, we used a simulated annealing rou- 
tine for the parameter search (Goffe et al. 1994), and an 
algorithm based on a combination of minimum square 
distance (in log L u ) between model and data, with the 
added constraint of reasonable agreement with the X-ray 


spectral index where available, enforced by increasing the 
score of the fit by a function of the angle between the 
model and observed spectral slope. 

We note that while simultaneous multi- wavelength ob- 
servations give the best ‘snapshot’ view of an entire SED, 
such data is usually still biased by being limited to a 
single epoch, as the variability of these objects is well- 
known. The compilation of all available photometry, gen- 
erally collected over a long time-frame, sampling many 
states of the SED, allows us to fit the model to the av- 
erage SED of the source which may be more appropriate 
in a population context, in which extreme source states 
might confuse results. 

2.4. The TEX Sample of Blazars 

After fitting with the SED model, each source was sub- 
sequently examined by eye, in order to identify ambigu- 
ous cases where the formal best-fit model did not clearly 
represent the only reasonable SED fit. In many cases 
these removed sources had few data points in the critical 
region from IR to UV, lacked X-ray spectral indices, or 
had significant host galaxy contamination in the optical. 
A comparison of an acceptable and a rejected source with 
their best SED fits is shown in Figure 3. A total of 653 
sources remained after this examination, in which we feel 
confident in the determination of the synchrotron peak 
to within one order of magnitude. (The rejected sources 
are discussed along with the other under-sampled blazars 
in §4.2.2). In total, 216 sources have robust estimates 
of ^ P eak» L peak 5 and L300; this set comprises the TEX 
(trusted extended) sample. Of these, 129 have extended 
radio flux from spectral decomposition, 87 were added 
using only VLA map data, and 90 have both a spec- 
tral and VLA estimate. For the remaining 437 blazars 
comprising the UEX (unknown extended) sample, upper 
limits on L 30 o (and subsequently Lki n ) were found by ex- 
trapolating from the luminosity at the lowest available 
frequency using a = 1.2. 

In Table 2 we list the sources which have a con- 
firmed synchrotron peak and extended power measure- 
ment (TEX sample). In columns 1 — 5 we give the 
source name, position, optical type, and redshift. Es- 
timated redshifts are noted by an asterisk in column 5. 
In columns 6 and 7 we give the best-fit synchrotron peak 
frequency and luminosity, respectively. In columns 8 and 
9 we give L 30 o and the converted value of Lkin, respec- 
tively. In column 10 we note the method for the extended 
radio flux estimate: VLA, spectral decomposition (SD), 
or both. We note those sources (19 total) in which the 
spectral result is used instead of a VLA measurement by 
a dagger in that column. In column 11 we give the ref- 
erences for VLA data. All values are in the rest-frame of 
the source. 

Because of the somewhat subjective requirement of an 
unambiguous synchrotron peak, the TEX sample does 
not contain all the sources in original list, and cannot be 
considered a complete sample. We focus most of our dis- 
cussion on qualitative analysis and positive findings for 
our well-characterized sample. Relaxing the criteria for 
inclusion in the sample risks inaccurate measurements of 
the three principle quantities (i/ pea k, L pea k> and Lkin) and 
could lead to incorrect conclusions. See §4.2 for a discus- 
sion of the selection effects of our master sample, the 
nature of the rejected sources, and the possible impact 
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TABLE 2 

TEX Combined Blazar Source List 


Name 

(i) 

R.A. 

(J2000) 

(2) 

Decl. 

(J2000) 

(3) 

Opt.* 

(4) 

z 

(5) 

log "peak 
(Hz) 
(6) 

log Lp ea k 
(ergs s -1 ) 

(7) 

log L 300 
(ergs s -1 ) 
(8) 

log L kin 
(ergs s^ 1 ) 

(9) 

Method 

(10) 

Ref. 

(ii) 

QSO 30003-06 

00 06 13.89 

-06 23 35.33 

B 

0.347 

13.16 

46.22 

41.68 

44.21 

both 

fjh 

QSO 30016+73 

00 19 45.73 

+73 27 30.07 

Q 

1.781 

13.26 

47.65 

42.61 

44.74 

both 

f 

QSO B0048-09 

00 50 41.31 

-09 29 05.21 

B 

0.200 

13.43 

44,97 

41.20 

43.93 

both 

b,d,f 

4C 01 02 

01 08 38.77 

+01 35 00.27 

Q 

2.099 

13.00 

47.71 

44.57 

45.85 

both 

c,f,h 

QSO 30109+224 

01 12 05.82 

+22 44 38.78 

B 

0.265 

13.87 

45.57 

40.33 

43.44 

both 

b,f,h 

QSO 30118-272 

01 20 31.66 

-27 01 24.65 

B 

0.557 

14.04 

46.13 

42.38 

44.61 

VLA 

d 

QSO Q0133+47 

01 36 58.58 

+47 51 29.30 

Q 

0.859 

13.34 

46.63 

41.93 

44.35 

VLA 

f,h 

QSO B0138 -097 

01 41 25.83 

-09 28 43.67 

B 

0.733 

13.49 

46.22 

42.41 

44.62 

VLA 

d 

4C 15.05 

02 04 50.41 

+15 14 11.04 

Q 

0.405 

12.39 

45.15 

41.39 

44.04 

VLA 

h 

QSO B0212+73 

02 17 30.82 

+73 49 32.55 

B 

2.367 

13.38 

47.20 

42.30 

44.56 

VLA 

f,h 

QSO B0215+015 

02 17 48.95 

+01 44 49.69 

Q 

1.715 

13.72 

47.71 

43.52 

45.26 

VLA 

b,f,h 

QSO B0218+357 

02 21 05.48 

+35 56 13.70 

Q 

0.960 

12.86 

45.91 

43.09 

45.01 

SD 


3C 66A 

02 22 39.61 

+43 02 07.80 

B 

0.444 

13.98 

46.34 

42.85 

44.87 

both 

b,e 


Note. — Table 2 is published in its entirety in the electronic edition of the Astrophysical Journal A portion is shown here for guidance 
regarding its form and content. Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 

Redshift estimated from optical magnitude 
i Spectral decomposition estimate of extended flux used instead of VLA estimate 
* Optical Spectral Type. B = BL Lac, Q — FSRQ. 
b Antonucci & Ulvestad (1985) 
c Murphy et al. (1993) 

' Cassaro et al. (1999) 
e Landt & Bignall (2008) 
f Kharb et al. (2010) 
g Condon et al. (2002) 
h Cooper et al. (2007) 

1 Caccianiga & March a (2004) 

on our conclusions. 

2.5. Radio Galaxies 

As the jets isolated in radio galaxies are understood 
to be the less-beamed counterparts to blazars, one 
way to understand the relationship of the synchrotron 
peak and beaming in the envelope is to see where the 
synchrotron peak falls for these jets. In the context of 
the presumed blazar sequence, these ‘hidden blazars’ 
should indicate whether the sequence is preserved in 
misaligned jets. We take advantage of recent work to 
isolate the radio, IR, optical, and X-ray core emission 
from the jet in radio galaxies using high-resolution 
mapping to place as many radio galaxies as possible 
on the i/peak - Lpeak plane (Chiaberge et al. 1999; 

Capet ti et al. 2000; Chiaberge et al. 2002; Capetti et al. 

2002; Hardcastle et al. 2003; Trussoni et al. 2003; 

Balmaverde et al. 2006; Massaro et al. 2010; 

Buttiglione et al. 2011). We examine the 45 radio 
galaxies with extended radio flux measurements, and at 
least one radio, optical, and X-ray measurement of the 
nuclear emission available in the literature. The L300 for 
all radio galaxies was found from a parabolic fit 2 to the 
steep, low-frequency radio emission and subsequently 
converted to L^i a as in Equation 1. In all cases where 
available, VLA map-derived estimates at 1.4 GHz were 
in agreement with the spectral fits. 

2 For fitting the steep component in radio galaxies, a log- 

parabolic model is more appropriate than a power-law due to the 
much greater range of the spectrum visible, which usually has some 
curvature in vFw 


Fourteen radio galaxies are nearby, well-known sources, 
with core measurements at multiple frequencies which al- 
low fitting of the synchrotron SED with the same param- 
eterization as for the blazars. This set includes seven of 
the Fermi-detected misaligned AGN (Abdo et al. 2010b). 

For the remaining 31 with sparser SED sampling, the 
synchrotron ^ pea k was estimated using a non- parametric 
likelihood estimator. The joint distribution of Rce> 
^peai.i and the SED colors 3 a ro , a ox was calculated for 
all the sources in our TEX, UEX, and fitted radio galaxy 
samples. From this, a conditional, one-dimensional den- 
sity on i/peak can be found by supplying the other three 
observables. The value of the peak was taken as the 
maximum of this distribution, and we also report the 
20% and 80% quartiles as the sampling error. 

It is straightforward to estimate the peak luminosity 
from the radio luminosity at 5 GHz, with which it has a 
linear correlation (correlation coefficient r = 0.85). The 
linear fit (with l-cr error) is 

log Lpeak =0.61 (±0.01) (log LgQHz — 43) 

± 45.68 ± (0.02) [ergs s -1 ] 

( 2 ) 

The 90% confidence interval, assuming a normal er- 
ror distribution, corresponds to an uncertainty of A (log 
Lpeak) — 0.8. 

The most likely peak frequency found by our statisti- 
cal model for all radio galaxies is at or below 10 14 Hz, 
with the exception of four FR II sources with moder- 

3 Spectral indices are defined at 1.4 GHz, 5000 A, and 1 keV. 
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TABLE 3 
Radio Galaxies 


Name 


R.A. 

Decl. 

Morph.* 

z 

log r'peak 

log Up ea k 

log L 300 

Uk.in 


(J2000) 

(J2000) 



(Hz) 

(ergs s _1 ) 

(ergs s _l ) 

(ergs 8“ 

(i) 


(2) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

3C 380* 

18 

29 

31.77 

+48 

44 

46.17 

2 

0.69 

12.95 

45.68 

44.06 

45.56 

3C 17 

00 

38 

20.52 

-02 

07 

40.8 

2 

0.220 

13.051°;?° 

45.35 

42.82 

44.86 

3C 18 

00 

40 

50.35 

+10 

03 

23.14 

2 

0.188 


44.79 

42.55 

44.70 

3C29 

00 

57 

34.92 

-01 : 

23 : 

27.89 

1 

0.045 

12.941®;!? 

44.01 

41.34 

44.02 

NGC 315 

00 

57 

48.88 

+30 

21 

08.81 

1 

0.016 

13.961®;!! 

43.94 

39.78 

43.13 

3C31 
3C 33.1 

01 

01 

07 

09 

24.93 

44.26 

+32 

+73 

24 

11 

45.20 

57.17 

1 

2 

0.017 

0.181 

13.091®!? 

12.991°;!! 

43.53 

44.30 

40.79 

42.42 

43.70 

44.63 

3C 66B* 

02 

23 

11.41 

+42 

59 

31.38 

1 

0.02 

13.36 

42.09 

40.73 

43.67 

3C 78* 

03 

08 

26.22 

+04 

06 

39.30 

1 

0.03 

12.99 

44.02 

40.94 

43.79 

3C 83. IB 

03 

18 

15.73 

+41 

61 

27.38 

1 

0.025 

13.09ig;|l 

43.40 

40.93 

43.78 

3C 84* 

03 

19 

48.16 

+41 

30 

42.10 

1 

0.02 

13.46 

43.09 

40.16 

43.34 

3C 111* 

04 

18 

21.27 

+38 

01 

35.80 

2 

0.05 

13.28 

43.29 

41.41 

44.05 

3C 133 
3C 135 
3C 165 

05 

05 

06 

02 

14 

43 

58.50 

08.34 

06.70 

+25 

+00 

+23 

. 16 
56 
i 19 

24.0 
31.63 
i 00.3 

2 

2 

2 

0.278 

0.127 

0.296 

13.04!°;|| 

13 -051®;| 

13.04!°?? 

45.08 

43.79 

44.33 

42.98 

42.04 

42.80 

44.95 

44.41 

44.84 

3C 171 

06 

55 

14.87 

+54 

08 

59.4 

2 

0.238 

13.041®;?? 

43.89 

42.73 

44.81 

3C 189* 

07 

58 

28.10 

+37 

47 

11.80 

1 

0.04 

13.44 

42.34 

40.88 

43.75 

3C 264 

11 

45 

05.00 

+19 

36 

22.74 

1 

0.022 

12.60!®!? 

43.77 

40.72 

43.66 

3C 270 

12 

19 

23.21 

+05 

49 

29.69 

1 

0.007 

1316i ils 

43.33 

40.09 

43.30 

NGC 4278 

12 

20 

06.82 

+29 

16 

50.71 

1 

0.002 

12.961®;?® 

42.55 

37.29 

41.70 

M 84* 

12 

25 

03.74 

+12 

53 

13.13 

1 

0.00 

12.98 

41.19 

39.04 

42.70 

M 87* 

12 

30 

49.42 

+12 

23 

28.04 

1 

0.00 

13.01 

41.58 

40.95 

43.79 

Cen A* 

13 

25 

27.61 

-43 

01 

08.80 

1 

0.00 

12.72 

41.47 

40.21 

43.37 

3C 287.1 

13 

32 

53.27 

+02 

00 

45.79 

2 

0.216 

13.05 +®S 

45.21 

42.38 

44.61 

3C 296* 

14 

16 

52.92 

+10 

48 

26.43 

1 

0.02 

12.99 

41.24 

40.57 

43.57 

3C300 

14 

23 

01.00 

+19 

i 35 

17.0 

2 

0.272 

12.731®;?° 

44.32 

42.85 

44.87 

3C 317* 

15 

16 

44.48 

+07 

01 

18.00 

1 

0.03 

12.83 

42.95 

41.34 

44.01 

3C321 
B2 1553+24 

15 

15 

31 

56 

43.47 

03.87 

+24 

+24 

04 

26 

18.93 

52.70 

2 

1 

0.096 

0.043 

13.271®;®° 

13.1lti:n 

44.19 

43.79 

41.86 

40.48 

44.31 

43.52 

3C 338* 

16 

28 

38.27 

+39 

33 

04.97 

1 

0.03 

13.12 

41.44 

41.18 

43.92 

NGC 6251* 

16 

32 

31.96 

+82 

32 

16.39 

1 

0.02 

12.79 

43.25 

40.68 

43.64 

3C 346* 

16 

43 

48.62 

+17 

15 

48.97 

1 

0.16 

13.46 

43.49 

42.25 

44.53 

Her A 

16 

51 

08.14 

+04 

59 

33.32 

2 

0.155 

13.041®;?? 

44.08 

43.60 

45.30 

3C 349 

16 

59 

32.00 

+47 

02 

16.0 

2 

0.205 

12.94:°;?? 

44.39 

42.49 

44.67 

3C 388 

18 

44 

02.41 

+45 

33 

29.81 

2 

0.091 

13.03!®;*! 

44.32 

42.00 

44.39 

3C 442A 
3C 449 

22 

22 

14 

31 

45.00 

20.62 

+13 

+39 

50 

21 

47.66 

29.81 

1 

1 

0.026 

0.017 

12.9812;?? 
!3.18l ?;| 

42.74 

43.27 

40.77 

40.17 

43.69 

43.35 

B2 2236+35 

22 

38 

29,46 

+35 

19 

40.01 

1 

0.027 

13.04+®;?? 

43.03 

39.62 

43.03 

3C 452 

22 

45 

48.78 

+39 

41 

15.36 

2 

0.081 

13.02!®;?? 

44.48 

43.34 

45.15 

3C 460 

23 

21 

28.50 

+23 

46 

i 48.0 

2 

0.268 

13.04!®;?* 

44.52 

42.52 

44.69 

3C 465 

23 

38 

29.39 

+27 

01 

53.52 

1 

0.031 

13.071®;?? 

44.07 

41.12 

43.89 


Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcmin- 
utes, and arcseconds. 

* Synchrotron peak estimated from blazar model fit to core fluxes 

* Morphological Type. 1 = FR I, 2 = FR II. 


ate jet powers (Lki n ~ 10 44,5 ergs s -1 ) which display 
extreme X-ray luminosities. These sources include the 
broad-lined radio galaxies (BLRG) 3C 390.3 (which has 
been previously reported to have an LBL-like SED by 
Ghosh & Soundararajaperumal 1995) and 3C 332, as 
well as two lesser-studied narrow line radio galaxies, 3C 
197.1 and 3C 184.1. The two BLRG have been de- 
tected with Swift and BeppoSAX in the hard X-ray 
(Grandi et al. 2006; Cusumano et al. 2010). All four 
sources were excluded from the final sample, as the X-ray 


is likely either thermal or from a strong IC component, 
and without knowledge of the synchrotron X-ray emis- 
sion, any peak estimate would be practically arbitrary. 

The final list of 41 radio galaxies is given in Table 3. 
Column descriptions are as in Table 2. 

3. RESULTS 

3.1. The Blazar Envelope 

The TEX sample of 216 blazars, and the 41 radio 
galaxies are shown in the plane of L pea k versus z^ pe ak in 
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Fig. 4. — The blazar sequence, which originally showed an anti- 
correlation between synchrotron L pea k and i/peak has been ex- 
panded into an ‘envelope’ with the addition of new observations 
and radio galaxies. BL Lacs are shown as filled circles, FSRQ as 
filled triangles, and radio galaxies as squares (FR I) and inverted 
triangles (FR II). Color indicates the jet kinetic power (L^in in ergs 
s -1 ), as estimated from extended radio flux measurements. Track 
(A) shows the path of a synchrotron peak for a single-component 
jet, and (B) for a decelerating jet of the type hypothesized to ex- 
ist in FR I sources. The fully-aligned limit for each (0 degrees) is 
shown as marked, with the arrow direction indicating the move- 
ment of the model source as it is misaligned. 

Figure 4. Blazars are represented as circles (BL Lacs) 
and triangles (FSRQ). The power of the source (Lion) is 
shown by color in bins of one decade. Examining the 
location of blazars, there is an obvious lack of sources 
with both high jet powers and synchrotron peak fre- 
quencies above 10 15 Hz, though a few of the ISP/HSP 
sources have quasar-like spectra. The population at low 
peak frequencies, alternately, is highly mixed in terms of 
kinetic power (with many low-power sources coincident 
with high- power sources), spectral type (with both BL 
Lac and quasar spectra present) and the spread of peak 
luminosities, which ranges from 10 44 to 10 47 * 5 ergs s -1 . 

The 41 radio galaxies are shown in Figure 4 as squares 
(FR I) and inverted triangles (FR II). Those detected 
by Fermi (Abdo et al. 2010b) are circled in red. For 
our radio galaxy sample, we note that the jet power 
does not apparently exert much influence on their syn- 
chrotron peak frequency, and it appears that all radio 
galaxies in our sample exhibit jets with peaks below the 
optical-UV range in frequency. Further, no evidence is 
found for any conservation of a ‘blazar sequence’ within 
the radio galaxy jets shown here. If, as has been sug- 
gested, different emission regions in the less powerful jets 
dominate at different orientation angles (Chiaberge et al. 
2000; Georganopoulos &; Kazanas 2003; Ghisellini et ah 
2005) it may simply be that the ‘sequence’ behavior (in 
which jet power regulates the synchrotron peak) pertains 
to the fast parts of the jet which are no longer visible in 
highly misaligned radio galaxies. 

It becomes interesting to ask the question of what 
type of orientation tracks are expected on the L pea k- 
^peak plane if we assume that powerful sources do not 


decelerate, while weak sources require velocity gradients. 
The former assumption is supported by the fact that 
the Lorentz factors required to model the blazar SED 
of these sources (T ~ 10 — 20) are similar with those 
deduced from following VLBI superluminal components 
(e.g. Jorstad et al. 2001), while the latter (as we dis- 
cussed in §1.2) is a result of our need to (i) explain the 
discrepancy between the high Lorentz factors required 
by the fast variability, pair production opacity, and SED 
modeling for the sub-pc scale jet and the lack of signifi- 
cantly superluminal motions at the VLBI scales and (ii) 
satisfy the unification between low power blazars and FR 
I radio galaxies. 

Two tracks are shown in Figure 4. These illustrative 
tracks do not intend to simulate any particular source, 
and are shown to simply demonstrate how sources that 
start with the jet pointing toward us will appear as the 
jet orientation angle increases. The first, labeled (A), 
follows the misalignment of a powerful source assuming 
a single jet Lorentz factor T = 15. This track follows 
the expected 1:4 ratio of decreasing (log) frequency to 
(log) luminosity of the peak as the source is decreasingly 
beamed (see, e.g., Lind & Blandford 1985). The second 
track, labeled (B), is for a weak source and shows the 
change in peak location for one possible model of a decel- 
erating jet, as described in Georganopoulos & Kazanas 
(2003). In this particular example, the Lorentz factor 
of the flow decreases continuously from T max = 15 to 
r m in = 3 in a distance of 5 x 10 17 cm, and the emitted 
spectrum is calculated along the jet and integrated spa- 
tially. The track exhibits the qualitative characteristics 
(e.g., more horizontal movement) not only consequent 
of the decelerating model, but of any model with veloc- 
ity gradients in which the most energetic electrons are 
injected in the fast part of the flow, such as the spine- 
sheath model of Ghisellini et al. (2005). Such models, as 
can be seen in Figure 4, exhibit a shallow decrease of 
Lpeak and a sharp decrease of i/ pea k because the opening 
angle of the beaming cone is smaller for higher frequency 
emission; as the angle to the line of sight increases, pro- 
gressively lower frequency emission is beamed away from 
the observer. Only for angles > l/r m i n , do we expect 
the beaming track to start resembling that of a constant 
Lorentz factor flow and L pea k to start dropping faster 
with decreasing y pea k- This model predicts that the low- 
power HSP blazars will follow a more horizontal path 
from the most beamed blazars to the misaligned radio 
galaxies. Interestingly, though each model originates in 
very different parts of the plane at low viewing angles, 
they coincide at higher orientation angles ( i.e when seen 
as radio galaxy core emission). 

3.2. Radio Core Dominance 

In order to investigate the effect of viewing angle in 
populating the blazar envelope, it is necessary to have a 
way to measure or approximate the degree of misalign- 
ment in individual sources. The most readily available 
measure is the radio core dominance Rce; typical values 
range from —3 (for the most off-axis radio galaxies) to 
~ 3 (highly beamed quasars). As the extended emission 
is non-variable and unaffected by beaming, the depen- 
dence of Rce on the orientation angle 0 is directly due 
to the increased Doppler boosting factor at low angles. 
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Fig. 5. — Core dominance (Rce) versus synchrotron peak lu- 
minosity for BL Lacs, FSRQ, and radio galaxies in our sample. 
Color indicates the kinetic jet power (Lkin in ergs s“ l ), as calcu- 
lated from extended radio flux measurements, as in Figure 4. The 
dotted line has a slope equal to 1.2, the OLS bisector slope found 
for the whole sample. This slope is consistent with the expected 
relation between L pea k and Lr, due to beaming effects. 

The Doppler factor 8 is given by 

6 = f (1 - PcosO) (3) 

where T is the Lorentz factor of the jet and P = v/c. The 
Doppler-boosted luminosity for any part of the jet SED 
is a function of 8 , 

L = Lq8 p (4) 

where L' is the unbeamed jet luminosity and p = 2 4- a 
to 3 + a depending on the geometry of the jet model, 
with a equal to the energy spectral index at the fre- 
quency of interest (see Urry & Padovani 1995). Rce is 
thus a monotonically decreasing function of 9, though 
the exact shape is dependent on T, becoming steeper as 
9 increases. This does create a difficulty in comparing 
Rce values directly for sources which may have veiy dif- 
ferent bulk T values (such as, perhaps, very high and low 
power sources). However, if we assume that sources of 
the same jet power are similar in T, then we can com- 
pare at least these directly, and in general a source with 
a lower Rce value will be more misaligned. 

In order to investigate further how the sources of simi- 
lar jet power are related in Figure 4, it would be helpful to 
add the core dominance. We thus examine how the core 
dominance changes with peak luminosity and frequency 
separately in Figures 5 and 6, still continuing to bin by 
the jet kinetic power as before. This choice is partly be- 
cause it is practically difficult to add a 4th parameter 
in Figure 4, and also because, as suggested in Figure 5, 
the core dominance may have a different maximum (i.e., 
normalization), depending on the jet power, such that 
high-power sources have a lower maximum Rce than that 
of the lowest power bin, making it difficult to compare 


sources at different powers directly. This can be under- 
stood if we make the assumption that unbeamed radio 
power scales from the overall jet power as = eLkin im- 
plying that jet radiative efficiency does not change much 
as the jet power increases. As we know Lkin— «Lf xt (tak- 
ing P ~ 0.6, see §2.2), this leads to the relation 

max Rce = (1 - 1/P) log L k i n -f p x log (2T) + c\ (5) 

where we have let <5 max =2r, c\ = log^ 1 /^), and (1-1/P) 
~ —2/3. Assuming the last term is roughly constant, 
the maximum value of Rce is thus expected to decrease 
with increasing L^n, even assuming reasonably different 
T values for high and low-power sources. For example, 
taking p=2, and T= 15 for a jet of Lkin = 10 46 predicts 
a max Rce about 1.25 lower than for a jet with Lkin — 
10 42 and T=3. 

It is also expected that sources will separate by jet 
kinetic power in Figure 5, as following from the above 
equations, we have for the expected relation between Rce 

and Lpeak^ 

Rce = log L p€ak + (l log L kin +c 2 (6) 

It is clear that binning by L k m will indeed induce the 
‘striping* in that figure through the second term, and 
that the slope in Figure 5 should be about 1.2 - 1.5 de- 
pending on the value of ctR (where we have assumed p = 
2 + a; see derivation in appendix B). However, scatter 
could be induced by the remaining term, C 2 = log^/c 1 /^) 
- (2-baj*)togLp eak /3 which is clearly determined by the 
relation between the unbeamed luminosity of a source 
at a given L k m- If this quantity is not roughly constant 
with respect to Lkin, Figure 5 would be reduced to a 
multi-color scatter plot and the striping would be lost 
or degraded. It is also interesting to note that, like in 
Figure 2, there is a fairly uniform, sharp transition from 
blazars to radio galaxies in Figure 5. These arguments 
again suggest that beaming has a dominant effect on the 
observed (core) luminosity, rather than being due to a 
range of unbeamed core luminosities (L^, Lp Cak ) for a 
given Lkin- 

A mild decrease in synchrotron peak frequency with 
increasing angle to the line of sight is expected in the 
simplest case of a jet dominated by a single Lorentz fac- 
tor. The change in peak frequency as core dominance 
decreases is not significant for our sample of sources with 
Lkin greater than 10 44 5 ergs s -1 , all of which have low 
synchrotron peak frequencies (these sources are therefore 
not shown in Figure 6). For low-power sources (less than 
10 44 ' 5 ergs s^ 1 ), however, there is a surprising pattern in 
the Rce — i'peak plane. As shown in Figure 6, there ap- 
pears to be a clustering of sources with low peak frequen- 
cies and yet high core dominance, where Rce reaches 
up to ~3, and another location at high peak frequencies 
where Rce reaches ~ 1—1.5. Either locus could conceiv- 
ably be connected via a misalignment track to the radio 
galaxies at lower left, but the source distribution sug- 
gests that there is not a smooth continuum between LSP 
and HSP sources. We expect a similar monotonic rela- 
tionship between Rce and i> pea k as was derived above 
for the peak luminosity, but this is clearly not consistent 
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Fig. 6. — Core dominance (Rce) versus the synchrotron peak fre- 
quency for all sources with L]^ below 10 44 - 5 (L 300 below ~ 10 42 ) 
ergs s _1 . While all low (kinetic) power radio galaxies have lower 
core dominance than blazars of similar power, there are apparently 
two locations of more aligned low-power blazars — one at low peak 
frequencies, with relatively higher Rce values (up to 3), and an- 
other at high peak frequencies, with slightly lower Rce values 
1—1.5). While nearly all sources at high peaks (> 10 15 Hz) are 
BL Lacs (i.e., HBLs), the low- peaking aligned sources include a 
mixture of FSRjQ and BL Lac spectral types. Color scale is the 
same as used in Figure 5. 

with what is shown in Figure 6, at least if sources at low 
power consist of a single population. 

In seme cases the high core dominance values for the 
LSP sources at low powers could be due to systematic un- 
derestimation of L 300 for these objects because of known 
problems with VLA measurements (incorrect modeling 
or use of an array configuration which might miss ex- 
tended flux). However this would require nearly every 
single source to be over-estimated by an order of magni- 
tude or more, and does not explain those LSP sources for 
which che steep spectrum is actually seen and matches 
the VLA estimate (such as BL Lac itself, noted in Figure 
6 ). 


4. DISCUSSION 
4.1. The Blazar Envelope 

Ghisellini et al. (1998) were able to reproduce the ob- 
served anti-correlations of the Fossati et al. (1998) blazar 
sequence using a simple one-zone leptonic jet model. 
They suggested that the increased radiative cooling of 
the electrons of high-power blazars lead to a particle en- 
ergy distribution with a break at lower energies (leading 
to the lower peak frequency). In support of the scenario 
of a continuous spectral sequence, our results in Figure 4 
clearly show a decreasing upper limit on L pca k as i/ pea k in- 
creases. That sources near this limit are highly aligned, 
and those falling below are progressively misaligned is 
supported by the strong correlation between core dom- 
inance and peak luminosity shown in Figure 5 (which 
we note is much harder to appreciate without separating 
sources according to their radio galaxy or blazar classi- 
fication and by their jet kinetic power). However, this 


simple scenario does not provide an explanation for the 
rather noticeable ‘L 1 shape to the source distribution in 
Figure 4 or the wide range of synchrotron peak frequen- 
cies for sources of similar low jet powers. While the for- 
mer might be the result of missing sources at intermedi- 
ate peak frequencies, the latter is a positive finding that 
cannot be explained away as a selection effect. 

4.1.1. Jet Structures in the Envelope 

Ghisellini & Tavecchio (2008) revisited the theoretical 
blazar sequence, linking the power of the jet and the SED 
to two fundamental parameters of the accretion disk — 
the black hole mass and the accretion rate. Following 
Narayan et al. (1997) they assumed that sources with an 
accretion power (Locc) lower than a critical value L cr = 
L EddV^cr are radiatively inefficient accretors with weak 
or absent broad line region, whiie sources with L acc > 
Lcr are characterized by efficiently radiating accretion 
disks with strong big blue bumps and powerful broad 
line emission (where L&dd is the Eddington luminosity 
and rher ~ 3 x 10” 3 — 10“ 2 a critical mass accretion rate 
in units of the Eddington accretion rate). 

As discussed in §3 we see in our results the manifesta- 
tion of two populations distinguished by typical jet struc- 
ture. The first population have non-decelerating jets 
characterized by a single Lorentz factor. These sources 
do not contain a high-peaking synchrotron component 
and the increase in peak frequency as the jet is aligned is 
modest, leading to a concentration of blazars with syn- 
chrotron peak frequencies ~ 10 14 Hz even when fully 
aligned. The sources in this population include all of 
the more powerful blazars (L^n > 10 44 ’ 5 ergs s -1 ) and 
many sources at lower powers. We include the latter be- 
cause of the high core dominance of many low-power, 
low-peaking blazars in Figure 6. The second population 
is only comprised of sources with Lki n below ~ 10 44 5 
ergs s” 1 and have complex, decelerating jets. These 
lower-power sources exhibit a high frequency-peaking 
synchrotron component when aligned. 

The non-decelerating, low-peaking (‘strong 1 jet) 
sources and decelerating, high-peaking (‘weak’ jet) 
sources form distinct populations which follow qualita- 
tively different de-beaming paths in L pea k — i/ pea k space, 
exampled by track (A) and (B) in Figure 4, respectively. 
The two population scenario explains the appearance of 
an L-shaped upper boundary as a natural consequence 
of the vertical and horizontal paths by sources through 
the L P eak — Z'peak plane. It also explains the appearance 
of Figure 6, as low power sources may be either strong 
or weak, producing considerable spread in possible i/ pea k 
values as Rce increases. 

It is tempting to associate those weak sources (having 
Lkin < 10 44 ' 5 erg s -1 ) that we argue are characterized 
by velocity profiles (track B in Figure 4) with sources 
having L acc < L^. If we assume that Lkin < L acc 
(e.g. Rawlings & Saunders 1991), this implies that the 
accretion power of these systems reaches at least up to 
Lacc ~ 10 44 5 erg s“\ which for thcr ~ 3x 10” 3 cor- 
responds to a black hole mass of ~ 10 9 M 0 , in agree- 
ment with black hole mass estimates (e.g., Barth et al. 
2002). Similarly, sources with L^n > 10 44,5 erg s -1 , 
which are all of the powerful strong jet type (track A 
in Figure 4), must have L aC c > L Ct > or equivalently 
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m > ihcr ~ 3 x 10“ 3 , assuming that black hole masses 
do not significantly exceed ~ 10 9 M 0 . 

In this framework, the low jet power sources (Lkin 
< 10 44 5 erg s^ 1 ) that yet seem to have a strong jet may 
be explained as sources with relatively smaller black hole 
masses, such that although their jets are low power, the 
mass accretion rate required to sustain them is m > m cr . 
Consider a jet with L^ n = 10 44 erg s" 1 and black hole 
mass 10 8 M 0 . From the above, it is clear that the ac- 
cretion rate (in units of 10” 3 ) m_ 3 = 0.7 Lki n ,44 Mg 1 
~ 7, exceeds the critical value placing this relatively low 
power jet in the family of ‘strong 5 jets. 

4.1.2. Optical Spectral Type and the Envelope 

While the original observational sequence of 
Fossati et al. (1998) implied a transition from quasars to 
line-less BL Lacs along a sequence from high to low bolo- 
metric luminosity, a firm distinction is not a necessary 
part of the theoretical scheme for a continuous blazar 
sequence. Ghisellini &; Tavecchio (2008) explain the 
possibility of blue FSRQ, formed when the dissipation 
radius for the jet exceeds the radius of the broad line 
region, as well as lower-power FSRQ formed when an 
efficient accretion forms around a relatively small black 
hole, (the explanation we proposed in §4.1.1 is along the 
same lines). However, the traditional divide between 
BL Lacs and FSRQ, at an equivalent width of 5 A may 
prove to be an essentially arbitrary one Ghisellini et al. 

a As discussed in Georganopoulos & Marscher 
, a line-less blazar may arise as a product of 
beamed emission out-shining the spectral lines, which 
would then emerge if the source were misaligned. A 
further complication arises when considering the spectral 
types of radio galaxies. High excitation radio galaxies 
(HERG) have prominent emission lines and presumably 
standard accretion disks accreting near the Eddington 
limit, while low excitation radio galaxies (LERG) have 
few or no lines and are likely radiatively inefficient. 
While most FR II are HERG and most FR I are LERG, 
there are known exceptions in both cases, including the 
nearby HERG FR I Cen A (Evans et al. 2004), other 
possible ‘FR I quasars’ (Blundell & Rawlings 2001; 
Heywood et al. 2007; Chiaberge et al. 2009), and an 
increasing population of LERG FR II at higher redshifts 
(up to z ~ 0.5, Garofalo et al. 2010). 

An explanation is given by Garofalo et al. (2010) for 
the presence of these mixed-type AGN based on black 
hole spin - namely, that central black holes with high 
retrograde spins and efficient accretion will produce pow- 
erful HERG FR IIs, which slowly spin down to prograde- 
spinning, low efficiency FR Is. However, if the transition 
between spin state occurs significantly before or after the 
transition in accretion efficiency, mixed-state objects will 
occur. Interestingly, the spin-based unification scenario 
may explain the lack of an obvious continuum of blazar 
types, as in this scenario the transition objects (near zero 
spin) have much lower jet powers than rapidly spinning 
FR I or FR IIs, thus in our typical luminosity-dependent 
plots there will be a ‘missing population’ of transition ob- 
jects, either showing up at much lower powers or missing 
due to selection effects. 

Ultimately, it is not clear what role the ionizing con- 
tinuum or line emitting regions has in influencing the jet 



Fig. 7. — The ratio log (L C ore/Lx) versus synchrotron peak 
frequency for all 653 sources with well-measured synchrotron 
peaks (equivalent a TX is shown on the right scale, a rx = log 
(L C ore/Lx)/7.68). This ratio is often used to estimate the location 
of the synchrotron peak frequency when detailed SED information 
is lacking. Landt & Bignall (2008) use log (Lcore/Lx) < 5 to indi- 
cate high synchrotron peak, however as shown, this ratio may only 
give a reliable upper limit to the frequency, and low log (L C ore/Lx) 
values do not correspond exclusively to high synchrotron peaks. 

type. Optical spectral type is therefore not a part of our 
suggested strong/weak classification, which relies instead 
on the source morphology, core-dominance, and location 
in the ^ P eak-L pea k plane. Thus, many of the low-peaking 
sources with high core dominance (t.e., strong jets) are 
BL Lacs, including BL Lac itself and all but one of the 
nine hot-spot BL Lacs found in the MOJAVE sample 
(Kharb et al. 2010). Similarly, several quasars appear at 
higher peak frequencies (i.e., as weak jets) in the UEX 
sample. 

Lastly, we comment on the possible presence of high- 
frequency-peaking FSRQ (HFSRQ) as claimed by sev- 
eral authors (Landt et al. 2006; Padovani et al. 2007; 
Landt et al. 2008). While the presence of quasars with 
low bolometric luminosities and high peaks (essentially 
having SEDs typical of HSP BL lacs) are consistent with 
expectations based on the updated theoretical sequence 
(Ghisellini & Tavecchio 2008), many of the objects dis- 
cussed do have higher jet kinetic powers typical of FR 
IIs. However, using low core radio to X-ray flux (log 
L C ore/log L x < 5) or high optical to X-ray ratio (a 0;r ) 
is not a reliable way to find sources with high i/ pea k. To 
illustrate this, in Figure 7 we plot log (L cor e/Lx), mea- 
sured at 5 GHz and 1 keV, against the known peak from 
our SED fits (based on much more extensive data than 
a two-point spectral index) for the 653 sources in our 
TEX-hUEX sample. While a few quasars exist in the 
HSP range > 10 16 Hz, these sources appear otherwise 
very similar to HSP BL Lacs in their SEDs. Further, 
many sources below the log (L core /Lx) < 5 limit clearly 
have low peak values, often with a hard X-ray spectrum 
indicating a non-synchrotron origin. So far no unambigu- 
ous evidence of any blazar with a high peak frequency 
and high luminosity (above the typical HSP zone) has 


The Blazar Envelope 


13 


TABLE 4 

Overview op Strong and Weak Jet Classes 



Strong 

Weak 

Definition 

More powerful, single-component jet (single char- 
acteristic T), remains relativistic on large scales 

Lower-power jet with faster deceleration, multiple 
emitting regions (multiple T) 

Optical Type 

mostly FSRQ 

mostly BL Lacs 

Synchrotron Peak 

All high-power LSP, some low-power LSP 

All ISP, HSP, some low-power LSP 

Typical L 300 

10 40 — 10 46 ergs s -1 

10 38 — 10 42 ergs s _1 

Typical L^in 

10 43 * 5 — 10 46 ergs s” 1 

10 42 — 10 44 - 5 ergs s -1 

Radio Morphology 

FR II; Collimated jets, hot spots, edge- brightened 

FR I; Less collimated, edge-darkened, no hot spots 

Accretion Type 

Efficient, likely standard disk 

Inefficient, highly sub- Eddington 


been found (see Maraschi et aL 2008; Landt et al. 2008). 

4.1.3. Intermediate and Low-peaking Blazars 

According to the proposed two-population scenario, all 
ISP and HSP sources belong to the weak jet parent sam- 
ple along with FR I radio galaxies. As blazars, they move 
from a population of HSP to ISP and finally LSP sources 
as they are misaligned. Contrary to the interpretation 
of a continuous sequence (Fossati et al. 1998), then, ISP 
blazars should have lower core dominance, lower overall 
luminosities, lower Doppler factors, and yet similar Lkm 
when compared with HSP sources. However, it is difficult 
to rule out a missing population of intermediate sources 
which may fill in the gap between low and high peaks in 
Figure 6 or in the center of Figure 4. These less popu- 
lated regions are not due to the flux-limited selections of 
the original master sample (see below). If these sources 
are in the master sample, the SED sampling would have 
to be much lower, and for such relatively bright objects, 
this seems unlikely. 

Weak jet sources at even greater angles to the line of 
sight should appear as LSP blazars, though they may not 
be recognizable as such (but rather as radio galaxies). 
Regardless, these should be distinguishable from strong 
jet LSP sources based on core dominance, morphology, 
extended radio power, and polarization signatures. How- 
ever without detailed VLBI and VLA studies to measure 
these attributes, confusion among the two LSP popula- 
tions is likely to occur. From our hypothesis that Rce 
must decrease with peak frequency for weak-jet sources, 
we expect most of the low-power LSP blazars shown in 
Figure 6 to belong to the parent sample of strong jet 
sources, along with FR II radio galaxies. This connec- 
tion is suggested by the equal or higher core dominance of 
most of these sources compared with the HSP blazars in 
the same Lki n range. It is also supported by the detailed 
observations we do have from VLBI studies of some of 
our objects, such as the nine LSP BL Lacs in the MO- 
JAVE sample with hotspots similar to those typically 
seen in high-powered quasars and associated with FR II 
morphology (Kharb et al. 2010). The MOJAVE blazars 
(almost entirely LSP sources) were also found to have an 
intrinsic, unbeamed luminosity function consistent with 
that of FR II radio galaxies, whether BL Lacs were in- 
cluded or not (Cara &: Lister 2008). 

We summarize the strong and weak jet characteristics 



Fig. 8 . — The blazar envelope for the UEX sample of 437 sources 
with only upper limits on Lki a . Upper limits were found by extrap- 
olating from the lowest frequency observation of core flux using a 
spectral index a = 1.2. The same color scale is used as in Figure 
4. We clearly see a similar ‘L’ shape as is suggested in the figure 
for TEX sources. Also shown are the three zones discussed in §4.2. 

as suggested by the current work in Table 4. 

4.2. Examining the Sample Bias 

Our effort is concerned with those sources which can be 
placed on the z/ pea k - L pea k plane with confidence, which 
is primarily a function of SED sampling. It is important 
to examine both the selection of the entire sample and 
the nature of the discarded objects for evidence of any 
bias which might produce the current findings. 

One possible selection bias is a ‘windowing’ effect due 
to the low frequency radio measurements required to iso- 
late the core and extended components. On one side, 
very core-dominated sources are likely to be missed when 
using only the spectral decomposition method, as these 
might not exhibit a steep spectrum except below very low 
frequencies (10 - 100 MHz) which are not frequently ob- 
served (and thus fail the two-component test). On the 
other side, a very de-beamed source may not qualify as 
a flat-spectrum source if the i / CTQ8 s is much above typ- 
ical sample definition frequencies of a few GHz. The 
combined effect could in principle lead to the selection 
of sources for which the core and extended fluxes are 
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by constraint rather similar. The inclusion of VLA map 
data, which allows us to find extended flux densities even 
in very core-dominated sources, should mitigate the first 
issue. As to the second, as we show in Figure 2, most 
blazars are core dominated above the crucial window of a 
few GHz, and we detect a few above this limit. Between 
the blazars and radio galaxies in this figure there is no 
large gap indicating a missing population, and we do not 
find any reasons to suspect that the windowing effect is 
a problem for the results of this paper. 

4.2.1. Upper Limits 

The over 400 sources in the UEX sample are shown 
in Figure 8, with the color corresponding to the upper 
limit of the jet kinetic power. The overall shape of the 
envelope appears qualitatively similar to Figure 4. 

We tested areas of the i+ e ak -L pea k plane using simula- 
tions to check how the sensitivity to sources in the plane 
is affected by our overall combined sample selection (ie., 
the combined selection of the samples in Table 1). We 
moved ~100 template SEDs drawn from our database to 
various peak locations in each zone shown in Figure 8 
and simulated cuts using the flux limits of our individual 
samples. In zones 1 and 2, we did not find any difficulty 
detecting the simulated sources out to a reasonable red- 
shift (z ~ 2). In zone 3, however, our sample would likely 
only include sources out to z ~ 0.2, which is not much less 
than the typical redshift of the HSP sources just above 
this region (z ~ 0.25). We conclude that the particular 
lack of ISP sources with L peak intermediate between LSP 
and HSP sources is not due to the initial selection of our 
sample. For at least zones 1 and 2, which would contain 
very bright sources, it is unlikely that a large number 
of sources would exist there, but lack the SED sampling 
necessary to make the more subjective cuts leading to 
our final sample. 

4.2.2. Remaining Sample 

We now examine the sources which we were unable to 
include due to inadequate SED sampling (typically fewer 
than ~ 10 SED data points). Examination of the colors 
(radio-optical and optical-X-ray indices) for these sources 
indicates that they have values consistent with the well- 
characterized sample of ~ 600 TEX+UEX sources. Un- 
surprisingly, those with X-ray detections mostly popu- 
lated the ‘HSP zone’ of the a ro — a ox plane, indicating 
that a large population of the excluded sources are likely 
higher redshift HSPs and ISPs. 

The optical and X-ray luminosity distributions for the 
excluded sample as a whole are indistinguishable (by KS 
test) from the TEX+UEX sample. In all three bands 
(and in particular the radio), those sources detected in 
the X-rays tended to have much lower luminosities, and 
are therefore unlikely to populate the so-called ‘forbid- 
den zone’ of high synchrotron peak luminosity and fre- 
quency. Those not detected in X-rays had a distribution 
of a ro typical of LSP sources, but actually tended to have 
higher luminosities than our TEX+UEX sample. How- 
ever, non-detection in X-rays and higher redshifts in gen- 
eral probably explains their absence in the TEX+UEX 
sample. Overall, none of the sources excluded from the 
TEX+UEX sub-sample appear to be substantially dif- 
ferent from those in it; however, only additional multi- 
frequency observations of these sources can confirm this. 


5. SUMMARY AND CONCLUSIONS 

In this work we have analyzed the multi-frequency 
spectrum of several hundred blazars in order to accu- 
rately measure the synchrotron peak frequency and lu- 
minosity, and performed spectral analysis of the low- 
frequency emission in order to estimate the intrinsic jet 
power by a scaling of the 300 MHz extended radio power. 
The connection of blazars and radio galaxies in the same 
plane amounts to an extension of the original blazar 
sequence to a blazar envelope, in which highly aligned 
blazars form an upper limit in the i/ pe ak — L pea k plane, 
while radio galaxies populate the lower left. We find that 
the jet power does not absolutely determine the ^ pea k 
of a maximally- beamed source. Instead, there seems to 
a break-down at low kinetic jet power, such that these 
jets are capable of having either low or high synchrotron 
peaks when aligned with the line of sight. 

Our results suggest a shift away from the continuous 
blazar sequence interpretation, but confirmation waits 
for ever more data. We summarize the most important 
points here: 

• Among our sample of over 600 blazars with well- 
sampled blazar SEDs, none with high peak frequen- 
cies and high luminosities are found. This void may 
be explained by the original arguments advanced 
by Ghisellini et al. (1998), that higher external ra- 
diation fields in high-power sources lead to greater 
cooling and thus lower break/peak energies. 

• However, it is not clear that there is a continuum 
of SED types with jet power. The source power 
alone does not appear to determine the placement 
of a source along the aligned sequence. Instead, we 
suggest that radio-loud AGN can be divided into 
two populations: low-synchrotron-peaking sources 
without significant velocity gradients or deceler- 
ation (strong jets) and high-synchrotron-peaking 
sources which do have such structure (weak jets). 
The attributes of these samples are summarized in 
Table 4. 

• It appears that below a certain jet power, 

sources can either have strong or weak jets, sug- 
gesting another parameter besides kinetic power 
is responsible for the jet structure. As in 

Ghisellini & Tavecchio (2008), we suggest that this 
parameter is the accretion rate m in Eddington 
units with the two families separated at m = 
m cr ~ 3 x 10 -3 - 10 -2 , with sources having 
rh < m cr exhibiting radiatively inefficient accre- 
tion (Narayan et al. 1997). 

• A key prediction of the two-population scenario is 
that ISP sources primarily consist of misaligned 
HSP objects, as suggested by their slightly lower 
Lpeak values and similar Lki n to HSPs. Observa- 
tions of ISP sources with synchrotron peak values 
and core dominance values between those typical 
of the LSP and HSP sources (currently missing) 
would support a continuous sequence interpreta- 
tion instead. 

• Observations of highly core-dominated, low- 
peaking BL Lacs (such as BL Lac itself) , including 
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several BL Lacs exhibiting FR II-like morphology, 
leads us to speculate that these sources belong to 
the strong-jet parent population, along with FR 
II radio galaxies. Likewise, some ISP and HSP 
sources with quasar-like spectra in our sample be- 
long to the weak population. Thus the spectral 
types associated with strong and weak jets are 
•nixed. 


• The low peak frequencies seen in our FR I radio 
galaxies supports the interpretation of a horizon- 
tal movement in the i^ pe ak — L pea k plane for weak 
sources which supports a model with velocity pro- 
files in the jet and leads to our two-population 
scheme. Observations of high synchrotron peaks 
in radio galaxies would be very difficult to explain 
under this scenario. 
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APPENDIX 

A. PHENOMENOLOGICAL BLAZAR SED MODEL 

The model was adapted from Fossati et al. (1997) and used to fit the average blazar SED with seven parameters, 
six of which were free to vary and which were optimized using an SA routine. 

Parameters: 


L r = 32 - 46 
r = 1000 

x j = 1 - 4 

«rx = 3-9 

a x = -0.4 — 0.8 

Speak = 10-20 

0 = 0.6 - 2.0 

To build the SED, the following equations were combined, where x — log(i7). Vector quantities are noted with arrows 

for clarity, and luminosities are in log. 


a = 



-^peak — P (^peak Xj 9.698) + Lr H~ 


C damp - { 2 _ 10 £i^ 
( P (x - 9.698) 4- Lr 


■f^synch — \ 


W - 


peak 


ddamp 


when x < ^ pea k 
else 

when X < (tfpeak - Xj) 

else 


Lie =L r -k ix + 7.685 4- (1 - a x ) [x - 17.383) 


(Al) 

(A2) 

(A3) 


(A4) 

(A5) 


Finally, we have the complete SED: 


L = log [l0^" ch + 10 £l °] 


(A6) 
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B. RADIO CORE DOMINANCE 

Equation 6 can be derived starting from the definition of log core dominance ratio Rce: 

Rce = log (-r^) (Bl) 

where L ex t and Lr are the luminosities of the extended and core components taken at the same frequency, and the 
beaming equations for the luminosity for the radio and synchrotron peak frequency are simply: 

Lr = Lr<5r" (B2) 

Lpeak = L' peak ^£ (B3) 

where 5 r and <5 pea k are the Doppler factors for the jet regions dominating the emission at those frequencies, and the 
exponents take values p = 24 -a or 3+ a, depending on the jet model, and with a the spectral indices in the band of 
interest. By taking advantage of the relationship Lkin = «Lf xt (see §3.2) we can rewrite it as: 

RcE=log .((^^) (B4) 

If we assume that = eLkm (see §3.2), and that the bulk Lorentz factors T for the radio and the peak emission are 
the same, such that 5r = <5 pe a k» using B3, we obtain: 

" ,<>g ( eL “- te) )-■<*(¥) i/# < b5) 

Letting pr — 2 4- an and p pea k = 2 + apeak = 3, we have equation 6: 

RcE - ( l°g Lpeak + (l - -p'j log Lkin + C2 (B6) 

c 2 = logleK 1 ^) - log Lpeak (B7) 
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